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^ . ABSTRACT 



' Context. The measurement of the scale of the baryon acoustic oscillations (BAO) in the galaxy power spectrum as a function of 

redshift is a promising method to constrain the equation-of-state parameter of the dark energy w. 

Aims. To measure the scale of the BAO precisely, a substantial volume of space must be surveyed. We test whether light-cone effects 
■ are important and whether the scaling relations used to compensate for an incorrect reference cosmology are in this case sufficiently 

accurate. We investigate the degeneracies in the cosmological parameters and the benefits of using the two-dimensional anisotropic 
power spectrum. Finally, we estimate the uncertainty with which w can be measured by proposed surveys at redshifts of about z = 3 
and z = 1, respectively. 

i-J^ ■ Methods. Our data is generated by cosmological N-body simulations of the standard ACDM scenario. We construct galaxy catalogs 

by "observing" the redshifts of different numbers of mock galaxies on a light cone at redshifts of about z = 3 and z = 1. From 
the "observed" redshifts, we calculate the distances, assuming a reference cosmology that depends on Wref. We do this for Wief = 
^ , -0.8,-1.0, and -1.2 holding the other cosmological parameters fixed. By fitting the corresponding (an)isotropic power spectra, we 

determine the apparent scale of the BAO and the corresponding w. 

Results. In the simulated survey we find that light-cone effects are small and that the simple scaling relations used to correct for 
the cosmological distortion work fairly well even for large survey volumes. The analysis of the two-dimensional anisotropic power 
spectra enables an independent determination to be made of the apparent scale of the BAO, perpendicular and parallel to the line of 
^SJ ■ sight. This is essential for two-parameter w-models, such as the redshift-dependent dark energy model w = wg + (I - a)w„. Using 

^ ' Planck priors for the matter and baryon density and AHq = 5% for the Hubble constant, we estimate that the BAO measurements of 

' future surveys around z = 3 and z = 1 will be able to constrain, independently of other cosmological probes, a constant w to ~ 12% 

l/^ I and ~ 11% (68% c.L), respectively. 
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1. Introduction ( 1200 ll) ). This has motivated consideration of more general dark 

. , ...... r , ■ r energy models that have a redshift-dependent EOS parameter 

S-^ There are two main observational indications of the existence of ^^^^ measurement of the parameter w can therefore help 

> an unknown energy component in the Universe, which is usually distinguish between not only a simple cosmological constant 

•"n referred to as dark energy (DE). First, observations or distant su- j »u j i j i t. » / n i u » »u 

S>y^ T r^. — — — T\ \, ^^J tT^ — H nrTTT;^ ^ and other dark energy models, but potentially also between these 

/S pernovae la (iRiess et alJfl998L iPerlmutter et al.lll999h favor an different models 
accelerating expansion of the Universe and therefore imply an 

energy component of negative pressure P. In particular, the pa- Several possible methods to constrain the EOS parame- 

rameter w of the equation of state (EOS) P = wp must obey ter w are summarized by the Dark Energy Task Force Report 

w < - 1 /3. Second, measurements of t he anisotropi cs in the cos- d Albrechtet al.||2006t ). The method that we consider here uses 

mic microwave background (CMB) dSpergel et al.ll2o5l l2007l) *e expansion history of the Universe. To measure this precisely, 

in combination with observations of the large-scale structure of require a standard candle or a standard ruler. For redshifts 

the Universe argue for a spatially flat Universe. Matter (bary- up to z ~ 1.8, supernovae la can be observed and calibrated 

onic and dark) contributes however less than 30% to the critical ^ be standard candles with which one can measure the lu- 

density. Hence, about 70% of the present-day energy density of mmo^^s distance Ddz). lEisenstein et al. | il998, 1999) proposed 

the Universe appears to be in an unknown form of energy. The *at the baryon acoustic oscillations (BAO) imprinted in the 

simplest way to account for this missing energy and the accel- S^iaxy power spectrum could be utihzed as a standai'd ruler, 

crating expansion is to introduce a cosmological constant A in The BAO have the same origin as the acoustic peaks in the 

Einstein's equations, which has a redshift-independent EOS pa- anisotropics of the CMB. Before recombination, baryons and 

rameter w = -1. So far all observations appear to agree with the photons were tighdy coupled. Gravitation and radiation pres- 

model of a cosmological constant. There is however at least one sure produced acoustic oscillations in this hot plasma; during 

theoretical drawback. The observed value of A, which is inter- *e expansion of the Universe, the plasma then cooled and fi- 

preted as vacuum energy, is highly inconsistent with cuiTcnt pre- "^lly nuclei and electrons recombined. The released photons 

dictions by particle physics, a discrepancy commonly ref erred to propagated through the expanding Universe and are observed 

as the cosmological constant problem (for a review, see cS^ by ourselves as the highly redshifted radiation of the CMB; 

in contrast, the baryons followed the clustering of the dark 

Senrf o^n«r re^Me.sfj to: C. Wagner, e-mail: cwagnerOaip . de matter and eventually collapsed to form galaxies. Since the 
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dark matter did not participate in the acoustic oscillations, the 
oscillatory feature in the galaxy power spectrum is far less 
pronounced than for the CMB photons (Peebles & Yu 1970: 
Sunvaev & Zeldovich 1970; Bond & Efstathiou 1984t|Holtzm3 
1989; .Hu & Sugivama.l996:.Eisenstein & Hu. 19981) . 



When the physical scale of the BAO has been calibrated 
using precise CMB measurements, it can be applied as a stan- 
dard ruler for measuring the angular diameter distance Z)/\(z) 
and Hubble parameter H(z). In contrast to supernovae la, the 
scale of the BAO is a more reliable standard ruler at high red- 
shifts. As the number of unperturbed peaks and troughs corre- 
sponding to the BAO in the power spectra increases, the wave- 
length of the BAO can be determined more accurately. On scales 
where structure growth is already nonlinear, the oscillations 
cannot be eas ily discer ned ("Spri ngel et al.ll2005t iGottlober et al.l 
[2006; Eisenst ein et al ]f2007b). Hence, with decreasing redshift 
the uncertainty in the observed scale of the BAO increases. To 
have good statistics for the first peaks a large volume has to 
be surveyed. The BAO have been detected in the present-day 
larges t galaxy redshift surveys (" Eisenstein et al.l2005l ; ICole et al.l 
l2005h . Data sets studying larger volumes and higher redshifts 
are however required to achieve tight constraints on dynami- 
cal DE models. Gal axy redshift surv eys at about z = 3, 1, or 



like HETDEX dHill et al. 120041). the NVFMO S BAO sur- 
('Glazebrook et al. 2005"; BassettetalJ HOGS'), and BOSS 



0.5, 

yey 

dSDSS-III Collaboration 2008), are designed for measuring the 
scale of the BAO and the EOS parameter w with a precision to a 
few percent and thereby constrain a variety of DE models. 

Many studies presented methods for extracting the scale 
of the BAO and estimated the accuracy of its measurement 
achievable b y future surveys. In these analyses, Monte Carlo 
simul a tions dBlake & Glazebrookl 120031: iGlazebrook & Blab ' 



20051; jBlake et al.' '2006'), Fisher matrix techniques ("Linde 
2OO3I; |Hu & Haiman 2003; Seo & Eisenstein 2003; Matsubara 
2004; Hutsi 2006b; Seo & Eisenstein 2007), N-body simula- 
tions (Angulo et al. 2005; White 2005- Seo & Eisenstein 200l 



iHuff et all l20 '07r iKoehler etaflT20()7t lArigulo et al.1 \200^. an d 
observational data dEisenstein et al.ll2005 : lHutsill2005Ll2006all3) 
were used. 

In this article, we attempt to include all important (physical 
and observational) effects for measuring the BAO and deriving 
constraints on the EOS parameter w. Our perspective is that of an 
observer, i.e. the starting point for the BAO measurement should 
be a galaxy catalog that provides celestial coordinates and red- 
shifts of the galaxies. Since we do not have in hand the type 
of real observations that we require, we have first to generate 
mock catalogs; we achieve this by completing "observations" 
on the data products of N-body simulations. Using these "obser- 
vations", we are able to study for the first time in detail the light- 
cone effect and the accuracy of the scaling relations used to com- 
pensate the cosmological distortion that arises by assuming an 
incorrect reference cosmology. A crucial point for all BAO mea- 
surements is the fitting method. We develop a method which uses 
only the oscillatory part of the power spectrum in a way that pro- 
duces unbiased and robust results. Further, we compare the re- 
sults of fitting the angle-averaged one-dimensional power spec- 
trum and anisotropic two-dimensional power spectrum. Finally, 
we predict the uncertainty with which proposed surveys will be 
able to measure the EOS parameter w by assuming two different 
w-models. 

The paper is structured as follows. In Sect. 2, we review how 
the EOS parameter w is measured using BAO. In Sect. 3, we 
explain how we generate the "observed" data from N-body sim- 
ulations, and in Sect. 4, we describe the power spectrum calcu- 



lation and our fitting method. In Sect. 5, we present our results, 
and finally we provide our conclusions in Sect. 6. 



2. Cosmology with baryon acoustic oscillations 

Knowledge of the true comoving scales, both parallel (ry) and 
transverse (r^) to the line of sight, of an observed physical prop- 
erty at a given redshift z, in the case of BAO statistical properties 
of large-scale structure, enables us to derive the Hubble parame- 
ter H(z) and angular diameter distance Da{z) from the measured 
quantities (redshift Az and angle A6) : 



ni = 



cAz 

W)' 

(I + z)Da{z)A6 . 



(1) 



In a flat Universe and for moderate redshifts at which the contri- 
bution of radiation can be neglected, the Hubble parameter de- 
pends on three cosmological parameters, the present-day Hubble 
parameter Hq, fraction of matter fin,, and the dark energy EOS 
parameter w(z), in the following manner 



H(z) = Hq + z)^' + (1 - n„) exp ^3 j| dzj . 

(2) 

The angular diameter distance Z)^(z) is a functional of the 
Hubble parameter and is given by 



Da(z) = 



1 +z 



dz 
H(z) 



(3) 



Assuming that we know apart from w all important cosmologi- 
cal parameters to significant precision from other observations, 
especially from CM B precision measurements such as WMAP 
dSpergel et al.l 120071) and direct measurements of the Hubble 
constant Hq by the HST Key project ( Freedman e t al. 2(301), one 
can constrain the EOS parameter w using this chain of equations. 
Since w appears in an integral over z in Eq. (|2|, observations at 
a single redshift cannot be used to measure the value of w at this 
redshift. A certain model for w has to be assumed before its value 
can be measured. In this paper, we assume a constant w and a 
simple redshift-dependent parameterization w - wq + (I - a)wa', 
this parameterization is used frequently in the literature (com- 
pare Dark Ene rgy Task For ce Report ( Albrecht et al. 2006)) and 
advocated by Lindeii d2007l) to be robust and applicable to a wide 
range of dark energy models. 

One problem of using the scale of the BAO as a standard 
ruler is that the BAO appear in a statistical quantity. We can- 
not measure Az of the BAO, but we can measure the redshifts 
of the galaxies and then, assuming a reference cosmology, re- 
construct their positions and derive their power spectrum. From 
this power spectrum, we can determine the apparent scale of 
the BAO and compare this with its true value. If they agree, 
we have used the correct cosmology. In principle, for every trial 
cosmology we have to recalculate the distances and recompute 
the power spectrum. A more efficient method is to scale ap- 
propriately the power spectrum derived for t he reference cos- 
mology using the following approximations (ISeo & EisensteinI 
2003;.Glazebrook & Blake 2005.) 



H'^'^iz) 



Da{z) ^ 



(4) 
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In this paper, we test if these approximations are still suffi- 
ciently accurate for very large surveys where z varies substan- 
tially within the survey. We achieve this by using three ref- 
erence cosmologies, which differ only in the EOS parameter 
Wref = -0.8, -1.0, -1.2. We use the relation (|4| to scale the 
derived power spectra and compare the results. For convenience, 
we define the following scaling factors parallel and transverse 
to the line of sight and an isotropic one for the one-dimensional 
angle -averaged power spectrum 

" W-Hz) ' ^ Da{z) ' (5) 

In Fig.[Tl we show the dependence of the scaling factors on a 
constant w for redshift z-'i (solid line) and z = 1 (dotted line). 
As a reference value, we chose Wief - -1.0. We also show the 
derivatives of the scaling factors with respect to w. The higher 
the value of the derivative, the more accurately we can constrain 
the constant w for a given uncertainty in the scaling factor. Since 
the derivatives decrease with decreasing w, the uncertainty in w 
will be, in general, higher towards lower w values. The derivative 
of /liso around w - -I at redshift z = 1 is ~ 1.5 times higher than 
at z = 3, that is although we can probably measure the scale of 
the BAO less accurately at redshift z = 1 than at z = 3 (due to 
nonlinear evolution), this does not need to be the case for the 
EOS parameter w. 




-1.5-1.4-1.3-1.2-1.1 -1 -0.9-0.8-0.7-0.6-0.5 -1.5-1.4-1.3-1.2-1.1 -1 -0.9-0.8-0.7-0.6-0.5 

W W 



Fig. 1. The scaling factors /i||_x and A\^a and their derivatives are shown 
as a function of w for the redshifts z = 3 (solid line) and z = 1 (dotted 
line). 



3. Simulation data 

To derive realistic samples of large galaxy surveys, we first com- 
pleted an N-body simulation in a large box (1.5 Gpc//i) and con- 
structed the corresponding dark matter distribution on a light 



cone by interpolating between about 20 snapshots. We calcu- 
lated the redshifts, which would be observed including the effect 
of peculiar velocities. Assuming a certain cosmology, we con- 
verted these redshifts into distances. Finally, we selected a cer- 
tain number of particles by applying a simple bias scheme and 
defining them as galaxies. In the following subsections, we de- 
scribe each step of this procedure. 

3.1. N-body simulation 

Our principal N-body simulation consists of 512^ dark mat- 
ter particles of a mass of 2 X lO'^M©//! in a (1.5Gpc//i) 
box. The initial power spe ctrum was produced by CMBfast 
(ISehak & Zaldarriagal Il996h . Starting from a glass distribu- 
tion, the particles were displaced according to second order 
Lagran gian perturbation theory by using the code of Sirl^ 
(l2005b . As cosmological parameters, we chose Qm - 0.27, 
Qa = 0.73, h = 0.7, Qb/i^ = 0.023, n, = 0.95, o-g = 0.8 
and w = -1.0. The simulation was performed with GADGET-2 
(Springel|20dl using a softening length of comoving 100 kpc//i. 
The starting redshift was z - 20. 

We also completed twelve 256^ dark matter particles simu- 
lations with the same cosmological parameters but different re- 
alizations of the initial conditions. We used these simulations to 
investigate systematic effects. Especially, we tested if our fitting 
procedure provides unbiased results. 

3.2. Ligiit-cone survey 

An observer at the present epoch fo (z(fo) - 0) identifies the 
galaxy distribution on his past light cone. He receives photons 
emitted at f < fo that have traversed the comoving distance 
X(t) = c J^'° a"'(f)df. To construct a light- cone survey, we follow 
an approach used bv lEvrard et al.l (|2002|) . We identify for each 
particle two consecutive snapshots between which it crosses the 
light cone and interpolate between them to find the position and 
velocity of that particle on the light cone. Expressed in formulas, 
the interpolated position is given by jc = jc,- + aAx where Xi is the 
position of the particle in snapshot / and Ax - Xi+i - Xi, and a is 
determined by requiring that \x\ - xiU + a At) with At being the 
time step between the two snapshots. After a Taylor expansion 
for the last term, we can solve for a 

a = , (6) 

2(xrAx+x(ti)Ax)' 

with Ax - x(fi) ~ > 0. The interpolated velocity is then 
given by V = (1 - a)vi + avi+i. 

For the galaxy sample at redshift z = 3, we place the cen- 
ter of the simulation box at redshift z = 3, i.e. 6.4 Gpc comov- 
ing distance away from the virtual observer. The orientation of 
the box is such that the line connecting the center of the box 
with the observer is parallel to the z-axis. The (1.5Gpc//i)^ box 
then extends from redshift 2.1 to 4.7, i.e. from 5.3 Gpc to 7.6 
Gpc in comoving distances. For this redshift range, we have 17 
snapshots at different times, i.e. of different expansion factors: 
a = 0.17, 0.18, .. .0.33. The box centered at redshiftz = 1 ranges 
from redshift 0.6 to 1.7. For this interval, we use 27 snapshots 
with a = 0.37, 0.38, . . .0.63 for the light cone construction. 

3.3. "Observation" 

Our virtual observer is at redshift z = and the simulation box 
is centered on redshift z = 3 and z = 1, respectively. To derive a 
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redshift for each particle, we compute the comoving distance x 
to the observer and solve the following equation numerically for 
the expansion factor a, 

r' da 

J a a^H(a) 

where H is the "true" Hubble parameter of the simulation 
H - H{) ^JQ.^/a^ + (1 - ^^m)■ The corresponding redshift is then 
Zh - 1 /a - 1 . In addition to this redshift caused by the Hubble ex- 
pansion of the Universe, there is also a Doppler redshift caused 
by the peculiar velocity of the particle. Hence, the total observed 
redshift is z = (1 + Zh)(1 + t^'rad/c) - 1, where Viad is the radial 
velocity of the particle. The accuracy of the simulation and the 
numerical algorithm of the integral equation above, implies an 
uncertainty of Az ~ 0.001. 

We next convert the observed redshifts to distances by as- 
suming different reference cosmologies, where the dark en- 
ergy EOS parameter w has the values w,ef = -0.8, -1.0, and 
-1.2. In this case, the Hubble parameter becomes H^'^^ia) — 

i/oV"m/a^ + (l-ilm)/a3<'+"-''. 



the M USYC collaboration (iGronwall et al.ll2007l: iGawiser et al] 
l2007h . 

For z - 1, we use the "weaker" bias sample, which is 
fairly consist ent with the expected bias of the target galaxie s 
of WFMOS ( Bassett et al.ll2005l:lThe KAOS Purple Boo]d2003h . 
The bias of the strongly biased sample at z = 1 is too high for 
these galaxies and we do not use this sample in the rest of the 
paper. 

3.5. Mock Catalogs 

For the survey centered on z = 3, we generate several mock 
catalogs of one million galaxies from the strongly and weakly 
biased light cones in redshift space using the entire 1.5Gpc/h 
box, which has a volume of lOGpc^. These numbers are at the 
upper limit of the current baseline of the HETDEX project. 

For the mock catalogs at z - 1, we use only the galaxies 
on the light cone around z = 1 in redshift space with a bias of 
b K 1.5. As a number of tracers, we choose one and two million 
galaxies in the entire box. 



3.4. Galaxy Bias 

Since the mass resolution of our simulation is too low for iden- 
tifying galaxy-sized friends-of-friends halos, we use a simple 
bias scheme to co mpile mock galaxy samples from th e dark 
matter distribution (ICole et al.lll998l:lYoshida et al.l200ll) . Using 
the initial density field, we define a probability function by 
P{v) cc (v-Vth)" if V is above the threshold Vth and zero otherwise. 
The dimensionless variable v is given by the density contrast nor- 
malized by its root-mean-square on the grid cr, i.e v{x) = S(x) / cr. 
The density contrast S(x) of each particle is computed in the fol- 
lowing way. First, we assign the particles with the cloud in cell 
(CIC) scheme to a 5 12^ grid dHocknev & Eastwoodll988D to ob- 
tain the density p on the grid points. Then, we calculate the den- 
sity contrast (5(j!:) = (p(jc)-p)/p on the grid points and interpolate 
it to the positions of the particles. 

In the next step, we Poisson sample the dark matter particles 
according to the probability function and track these "galaxies" 
throughout the snapshots and the constructed light cones. 

As parameters, we use (a = 0.3, Vth = 0.4) for a strongly 
biased sample and (a = 0.2, Vth - -0.6) for a more weakly bi- 
ased sample. We calculate the corresponding galaxy bias to be 
the square root of the ratio of the galaxy power spectrum to the 
dark matter power spectrum in real space: = Po-i\{k)l Pdmik). 
The bias of the mock galaxies decreases with decreasing redshift 
(Fig|2|l as expected in the model. We select primarily particles in 
high-density regions of the initial density field, which occupy 
less prominent structures at later times due to the further devel- 
opment of gravitational clustering. The bias has a mild scale de- 
pendence due to the schematic procedure of grid-based density 
estimation and the specific probability selection function. 

The parameters a and Vth were chosen such that the strong 
bias sample at redshift z = 3 is consist ent with the ex pected 
bias of the target galaxies of HETDEX (Hill et all 12004 ). The 
bias of this sample is very similar to the friends-of-friends halo 
sample with a halo mass higher than ~ 4 x 10'^M(T)//i ob tained 
from the Mare Nostrum simulation (iGottlober et al.l 12006) . As 
observations of the target galaxies of HE TDEX, namely Lyman- 
a e mitting galaxies (LAE), suggest (" Gronwall et al.l 120071: 
IGawiser et al. 2007), this expected value for the bias might be 
too optimistic. Therefore, we also use the weaker bias sample at 
z = 3, which matches the measured bias of LAE observed by 
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Fig. 2. Galaxy bias for the two different parameter sets at redshift z = 3 
and z= 1. To construct this plot, we used 15 million galaxies to reduce 
the shot noise. 



4. Analyzing the data 

We analyze the mock catalogs first by calculating the power 
spectrum and then by fitting the extracted BAO. In the follow- 
ing two subsections, we describe how we calculate the power 
spectrum and provide details about the fitting procedure. 

4.1. Power spectrum 

Since we use an FFT to complete the Fourier transformation, 
we first have to assign the particles to a regular grid. We con- 
sider a 1024^ regular grid and select the CIC scheme to com- 
plete the mass assignment. After converting the density field 
p(x) to the density contrast 6{x) - ipix) - p)lp and performing 
the FFT to evaluate the Fourier transform (FT) 6(k), we com- 
pute the raw isotropic (one-dimensional) and anisotropic (two- 
dimensional) power spectrum by averaging |<5(^)|^ over spher- 
ical shells PZ^ik) = (|(5(^)l')^j^^„ and rings P-;-(^||,^J = 
^|<5(^)P^ . , respectively. These raw power spectra are related 
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to the true power spectrum in the following way ( lJindl2005l) : 



piso./a„iso. ^ / 2 |^(^ ^ 2k^^mtP{k + 2k^ym) 

+ ^J]mk + 2k^ym)A , (8) 



meZ' 



shell/ring 



where W{k) is the FT of the mass assignment function, the 
Nyquist frequency, and n the number density of tracers for the 
density field. The first term in the equation above reflects the 
convolution of the density field with the mass assignment func- 
tion. The second term corresponds to the shot noise. The sum- 
mation over all integer vectors m takes care of the aliasing. In 
the case of the CIC scheme, W(k) takes the following form 



W(k) 



sin 

n- 



(9) 



To derive a good estimate of the true power spectrum from 
its raw form, we first subtract the shot noise and then follow an 
iterative method for the correction of deconvolution and alias- 
ing as proposed by Jing (2005). We subtract only for the galaxy 
samples shot noise, since we do not detect any shot noise for the 
dark matter densities. 

We estimate the error in the power spectrum by count- 
ing the number of independent modes used in the calcula- 
tion of Piaw, which corresponds approximately to cr™° - 



(P + Ns) and cr 



.shell 



Ak- 



■ (P + Ns) 



jFeldman et al.l [T9941) for the anisotropic and isotropic power 
spectrum, respectively. Here, Nk. denotes the thickness of the 
rings and shells we averaged over For the anisotropic power 
spectrum, there is additionally the parameter AyU which is the 
range in the values of the cosine of the angle between the wave 
vector and the line of sight. The shot noise is denoted by A^,. All 
the power spectra in this paper have independent bins with a bin 
width of A;t = 0.005 /i/Mpc. 




k [h/Mpcl 

Fig. 3. Power spectrum of the dark matter distribution at z = 3 (lower 
curves) and z = 1 (upper curves). The dashed lines are the linearly 
evolved initial power spectrum. In the case of redshift space, the linear 
power spectrum was multiplied in addition by 1-1-2/3 fi{z) +1/5 /?(z)^ ~ 
1.843(1.714) forz = 3(1). 



z=3 real space 
z=3 redshift space 

z=1 real space 
z^1 redshift space 
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Fig. 4. Fractional dilference of the light-cone power spectrum and 
the power spectrum at the corresponding snapshot Psnap- The horizontal 
lines show the predicted shifts from linear theory in real space (dashed) 
and redshift space (solid). 
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Fig. 5. Fractional difference of the galaxy light-cone power spectrum 
Pic and the galaxy power spectrum at the corresponding snapshot Psnap- 



In Fig. |3] we show the isotropic dark matter power spectra 
at redshift z - 3 and z - I obtained by assuming the correct 
cosmology with Wref - -1. If not stated otherwise, we always 
use as the reference cosmology the cosmology of the simula- 
tion, i.e. we assume that Wmf - -I - The linear redshift distortion 
(Kaiser 1987) amplifies the power by 1 + 2/3y6(z) + l/5/3{zf 
with /3{z) = f(z)/b{z), where f{z) is the linear growth rate and 
b(z) the bias parameter. In the case of dark matter we have 
b - I. For scales shown, the nonlinear evolution is still mild. 
The deviations in the real-space power spectra from the lin- 
early evolved initial power spectra (plotted as dashed lines) start 
around k ~ 0.25/!/Mpc and k ~ 0.15/!/Mpc for z = 3 and z = 1, 
respectively. The nonlinear redshift distortions ("finger of God" 
effects) in redshift space lead to a suppression of power on small 
scales, which almost balances the increase in power due to non- 
linear clustering at z = 1 (filled squares in Fig.[3]l. The BAO can 
be seen as tiny wiggles in the power spectrum. Their amplitudes 
are < 5% of P{k). 

We do not show the light-cone power spectra, since they 
lie almost exactly on the corresponding snapshot power spec- 
tra. The fractional differences in the power spectra derived from 
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the light cone and the corresponding snapshot are instead plotted 
in Fig.|4] On linear scales, the fractional differences are ~ 0.5% 
and ~ -2.5% at z - 3 and z = 1, respectively. The reason for 
these differences is that the value of the growth function at the 
mean redshift is not equal to the mean growth function. The 
stated numbers can be understood by averaging the square of 
the redshift-dependent growth function multiplied by an appro- 
priate geometrical factor over the survey box (Yamam oto et al.l 
119991) . Numerical calculations for our survey designs are shown 
as dashed (real space) and solid (redshift space) lines. The de- 
viations at larger k are due to nonlinear effects. The question of 
whether these light-cone effects alter the fitting of the BAO is 
addressed in Sect. |5] 

In Fig. |5] we compare the biased galaxy power spectra de- 
rived from the light cones with those evaluated for the corre- 
sponding snapshots. Qualitatively, we observe the same behavior 
as in the dark matter case. The larger scatter is due to higher shot 
noise and the means are shifted slightly, since, as happened be- 
fore for the growth function, the bias at the mean redshift is not 
equal to the bias averaged over the appropriate redshift range. 



4.2. Fitting procedure 

Our fitting method attempts to remove all the information apart 
from the BAO themselves. This is easily achieved by dividing 
the power spectrum by a smoothed version of itself. The advan- 
tage of this is that one does not need to model all the physi- 
cal processes that alter the power spectrum such as redshift dis- 
tortion, nonlinear evolution, and galaxy bias or uncertainties in 
cosmological parameters, such as erg, «s, and massive neutrinos, 
which affect the overall shape of the power spectrum but have 
little effect on the BAO. 

The extracted BAO of the reference power spectra were then 
fitted to the extracted BAO, allowing the scaling factors A^^ or 
/liso to vary. In an additional fitting attempt, our free parameter is 
instead the EOS parameter w and we derive the scaling factors 
/t|U or '^iso by applying Eq. dSjl. 

The fitting parameters are determined by Monte Carlo 
Markov chain (MC MC) techniques. We use the Metropo lis- 
Hastings algorithm (iMetropolis et al.lll953t lHastings|[T970h to 
build up the Markov chain. 

In the following, we describe each step of the fitting proce- 
dure. 



Smoothing 

There are previous proposals to use only the oscillatory 
part o f the power spectrum for r neasuring the scale of the 



BAO ("Blake & Glazebrookl 120031: IGlazebrook & Blakel 
Hutsi 2005, 2006a; Koehleret2]Bo07; Anguloetal.' 



1 


2005 




2008 



Li 

IPercival et al.ll2007 b). All of them divide the power spectrum by 
(or subtract from it) a non-oscillating fit. The methods used to 
derive the smooth fit range from deriving a semi-analytic zero- 
baryon reference power spectrum, fitting with a quadratic poly- 
nomial in log-log space, or using a cubic spline, to fitting with 
a non-oscillating phenomenological function. In this paper, we 
generate a smoothed power spectrum in an almost parameter- 
free way, by computing for each point the arithmetic mean in 
log space of its neighbors in a range of +0.03 /i/Mpc in k. In the 
two-dimensional case, we smooth radially in the direction of |^|. 
This smoothing length is the only parameter in our smoothing 
method, which, for all input power spectra, we select to have 



the same value of approximately equal to the wavelength of the 
BAO. 

By dividing the measured power spectrum by its smoothed 
version, we derive the purely oscillatory part of the power spec- 
trum Po.sci = P/Psmooth (see Fig.©. 

Fitting and Priors 

In pl ace of fitting the extracted BAO by a (modified) sine fu nc- 
tion dBlake & Glazebrook 2003: Glazebr ook & BlaSl2005h or 
using a periodogram (Hiitsi 2005) to measure the scale of the 
BAO, we compare the BAO with a range of different oscillatory 
reference power spectra PJ^'^^ produced from thousands of lin- 
ear power spectra generated with CMBfast, which differ in the 
cosmological parameters O^, fib, Hq, and «s. 

We prefer this method since the BAO are not exactly har- 
monic dEisenstein & Hulll998l;lKoehler et al.ll2007l) and the un- 
certainties in the aforementioned cosmological parameters can 
easily be included. As priors on these cosmological parame- 
ters we use th e predicted uncertainties for the Planck mission 
(iThe Planck Bluebook 2005.) . For the Hubble constant Ho, we 
combine the P lanck priors with con straints provided by large- 
scale structure (IPercival et al ."2007a') and include measurements 
from the HST Key project (jpreedman et al. 2001); 



Aw,n = 1.25% , AHo = 5% , 
Awb = 0.75% , An, = 0.5% 



(10) 



where cj^ = Q.^ and = h^- We assume throughout this 
article that the Universe is spatially flat. 

We determine the scaling factors /iy x or /liso by fitting a 



scaled P'^ljMh or P^M^^ok) to P,,,-,. 

Damping and Reconstruction 

Since nonli near structure growth di minishes the amplitudes of 
the wiggles dEisenstein et al.ll2007bl) . we can improve the fitting 
by adding a suppression factor W = exp(-crA;^) to our fitting 
function 

Posciik) = W [Pltii^i^ok) - l] + 1 . (11) 

For the two-dimensional redshift-space power spectrum, we 
have to introduce two suppression parameters o-y^, since the 
redshift space distortion increases suppression along the line of 
sight. The suppression factor then becomes W = exp(-cr||A;y - 

(r±kl). 

If the density field is known accurately, i.e. the shot noise is 
small and the galaxy bias is k nown, it is possible t o use recon- 
struction techniques (see e.g. iNaravanan & Croffl d 19991) for a 
comparison of methods) to undo, at least in part, the suppression 
of the wiggles (Eisenstein et al. 2007a). We applied the PIZA 
method (Croft & Gaztanaga 1997) to the dark matter distribu- 
tions. We attempted a similar technique for the mock catalogs 
but due to shot noise and galaxy bias our results were unsatis- 
factory. Further efforts are required to achieve this goal. In Fig. 
|6] we indicate the BAO extracted from the data by crosses. The 
dashed lines show the best-fit functions and the dotted lines de- 
pict the corresponding non-damped reference P^^^^.iik). We ob- 
serve that the suppression of the BAO is higher in redshift space 
and increases with time. Additionally, the results of the recon- 
structed density fields are plotted (data as filled circles and the 
best-fits as solid lines). The amplitude of the BAO could not be 
reestablished completely but to a significant fraction. 



Wagner et al.: Constraining DE via BAO in the (an)isotropic light-cone power spectrum 



7 



1.10 



1.05 



E 1.00 



0.95 



0.90 



T 1 1— 



z=3 real space 
reconstructed 




1.10 



1.05 



E 1.00 



0.95 



0.90 



T 1 1 — 



z=1 real space 
reconstructed 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
k [h/Mpc] 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
k [h/Mpc] 



1.10 



1.05 



E 1.00 



0.95 



0.90 



z=3 redshift space 
reconstructed 




1.10 



E 1.00 



T 1— 



z=1 redshift space 
reconstructed 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
k [h/Mpc] 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
k [h/Mpc] 



Fig. 6. The left upper panel shows the oscillatory part of the power spectrum at redshift z = 3 obtained from real space and its best-fit function 
including the suppression factor (dashed line). The dotted line corresponds to the oscillatory part of the initial power spectrum. The solid line is 
the best-fit of the reconstructed data. At the bottom left, we present the same in redshift space. On the right the corresponding plots for z = 1 are 
shown. 



Fitting Range 

The k range over which we use the data for our fitting plays 
an important role. In the fitting procedure we use a minimum 
wave number fcmin = 0.04/!/Mpc, since for the points with 
k < 0.03 h /Mpc the smoothing is not well defined, and the error 
due to sample variance is high. As the maximum wave number, 
we select k^^^ = 0.25 (0.30), for redshift z = 1 (3). Although 
nonlinear evolution has already started on these scales, the BAO 
are still clearly visible (see Fig.|6]l. 

5. Results 

The most important measurements in our analysis are the 
marginalized probability distribution functions (PDF) of the 
scaling factors and the EOS parameter w. To determine their val- 
ues, we marginalize the joint PDF produced via the MCMC tech- 
nique over all other fitting parameters. From these functions, we 
can derive the best-fit values and the accuracy of both the scaling 
factors and the EOS w. 

Fitting mettiod, real space versus redstiift space, and degen- 
eracies 

With the twelve 256"' simulations we assess the robustness of our 
fitting method and the presence of systematic effects. We deter- 



mine the best-fitting scaling factors, both parallel and transverse 
to the line of sight, for the dark matter power spectrum in real 
and redshift space for all twelve simulations at redshift z = 10, 
3, and 1. For all redshifts, we find that the mean value of the 
twelve simulations is clearly less than Icr (standard deviation) 
away from unity. Hence, our fitting method provides unbiased 
results within the margins of error (see Table[T]i. 



Table 1. Results for the scaling factors /ly and Aj_ obtained from the 
dozen 256^ simulations. 



Parameter 


Real Space 
mean sigma 


Redshift Space 
mean sigma 


Difference 
mean sigma 


z= 10 

-ill 


0.998 0.007 
0.999 0.005 


0.998 0.008 
0.999 0.005 


0.000 0.001 
0.000 0.001 


z = 3 
-ill 


0.998 0.009 
1.000 0.006 


0.999 0.011 
0.999 0.006 


-0.002 0.004 
0.001 0.002 


z= 1 

-ill 

-Ix 


0.997 0.010 
0.999 0.009 


0.999 0.014 
0.999 0.010 


-0.003 0.008 
0.000 0.003 



The differences between results for real and redshift space 
in each simulation do not show a systematic trend. The mean 
difference is in fact consistent with zero. The error in the scale 
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factor parallel to the line of sight A\\ is larger in redshift space, 
since the damping of the BAO is stronger in redshift space (see 
Fig©. 

To understand the degeneracies in the scaling factor for the 
cosmological parameters included in our method, we complete 
the fitting by keeping all cosmological parameters fixed to the 
value of the simulation apart from the parameter under consider- 
ation, for which we calculate the PDF of Also for different values 
of this parameter. The corresponding 95% confidence lines for 
the dark matter power spectrum at z = 3 (solid) and at z = 1 
(dotted) are shown in Fig. ITH] Since the sound horizon imprinted 
in the power spectrum is redshift independent, we expect that 
the derived degeneracies are also redshift independent. This is 
indeed the case for our fitting method. The only difference with 
respect to z is the larger uncertainty in Ai^o at lower redshifts due 
to the suppression of the BAO by nonlinear evolution. 

The interval of cosmological parameters plotted in Fig.|2]was 
chosen to be centered on the values of the simulation and in the 
range of +3%. In this way, we can immediately assess the degree 
of degeneracy. We find that if we alter one of the cosmological 
parameters by 1% the scaling factor changes by 0.22% for w^, 
0. 14% for (jjb, ~ 0% for Ho, and 0% for n,. These numbers agree 
w ell with those given by th e fitting formula for the sound horizon 
bv lEisenstein & Hul (Il998h : 



44.51n(9.83/wn,) 



1 -(- lOw, 



3/4 



Mpc. 



(12) 
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to z = 3 and the dotted line to z = 1 . In contrast to the previous 
case, these correlations are redshift dependent, for example at 
z = 1 the dependence of both the sound horizon and the Hubble 
parameter H(z) on w,n almost exactly cancel each other, whereas 
at z = 3 the effect originating in the Hubble parameter H{z) is 
stronger. Among the cosmological parameters the largest uncer- 
tainty in Wo originates in //q, in particular at lower redshifts. 




1.10 -1.05 -1.00 -0.95 -0.90 



Fig. 7. Dependence of the fitted A[so on the cosmological parameters 
Wm, cob, Hq, and at redshift z = 3 (solid) and z = 1 (dotted). 



Fig. 8. Same as in Fig.|7]but for wq. 



Light-cone effect 

We have already mentioned that the power spectra derived from 
the light cone output and corresponding snapshot at the mean 
redshift are almost on top of each other. The more evolved parts 
of the light-cone sample compensate almost precisely the less 
evolved parts, such that its power spectrum is almost identical 
to that at the mean redshift (see Fig. HJl. In Fig. |9]we show the 
corresponding PDFs of the scaling factors. The lines show the 
68% and 95% confidence levels obtained from the dark matter 
snapshot (dashed) and light-cone (dotted) power spectra at red- 
shift z = 3 (left) and z = 1 (right) in real (top) and redshift space 
(bottom). The differences in the error ellipses derived from the 
light-cone and snapshot data are overall small. Hence, light-cone 
effects in this survey volume and at these redshifts are unimpor- 
tant to our fitting. 



More interesting are the dependences of the equation of state 
w on the cosmological parameters, in particular the dependence 
on oJtn and //q, which enter the relation between the scaling fac- 
tors and w. For a constant w model (w - wq), we indicate the 
derived correlations in Fig. [8] where the solid line corresponds 

' Since we used only one realization (the 512^ simulation) to produce 
these plots, the fitting results were slightly off the input values of the 
simulation. For display purposes, we centered the lines. 



Reconstruction 

The error ellipses for the reconstructed samples are shown in 
Fig.|9]as solid lines. We observe that for data at z = 1 in redshift 
space the reconstruction of the BAO shrinks the error ellipses by 
a factor of three. For this reason, it would be desirable to develop 
a reconstruction method that can be applied to noisy and biased 
density fields. For surveys at redshift z = 3, this effect is less 
important. 
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Orientation of the error ellipse 

The errors in the scaling factors originate in two different 
sources. One source is the errors in the power spectrum; the 
other source is the uncertainties in the cosmological parameters 
and which produce an uncertainty in the sound horizon, 
i.e. the scale of the BAO. If the principal error could be attributed 
to the uncertainty in the physical sound horizon, the orientation 
of the ellipse would be approximately in the direction of the line 
defined by = For a low signal-to-noise ratio power spec- 
trum for which the scale of the BAO in the power spectrum can- 
not be determined accurately, the orientation of the error ellipse 

is instead in the direction of the line of ^A]_A\\ = const. 





0.96 0.98 1.00 1.02 1.04 
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Fig. 9. The marginahzed probability distribution functions (PDF) of the 
scaling factors at redshifts z = 3 (left) and z = 1 (right) in real space 
(top) and redshift space (bottom) derived from the light cone (dotted), 
snapshot (dashed), and "reconstructed" light cone (solid). 




0.15 0.2 
k [h/Mpc] 



Fig. 10. The oscillatory part of the real space power spectra at red- 
shift z = 3 derived for three different cosmologies and its best-fit func- 
tion. For display purposes, the power spectrum for Wref = -1.2(-0.8) is 
shifted -l-0.2(-0.2) along the y-axis. 



the mean redshift, which is not equal to the mean scaling factor 
The noise in the scaled power spectra is substantially smaller 
on large scales than the error, which is already implicit in the 
power spectrum due to cosmic variance (dashed line) and ad- 
ditional shot noise (solid line). To determine if this (additional) 
noise impairs the fitting of the BAO and with it the measure- 
ment of the equation of state w, we fit the corresponding data of 
the twelve 256-' simulations by assuming a constant w. We find 
that the best-fit value for wq varies among the different reference 
cosmologies for a single simulation by not more than 4%, and 
is typically ~ 2%. This should be compared with the standard 
deviation for a single wq measurement, which is about 3%. All 
PDFs scatter about the mean value of wo = -1, no systematic 
effects are noticeable. We conclude that the scaling relations do 
not introduce a noticeable bias or enlarge the error in w. For real 
data, a consistency check would be to recalculate the distances 
and corresponding power spectrum and redo the analysis using 
the measured w as the reference value for the assumed cosmol- 
ogy. 



Cosmological distortion 

We assess if the approximations to compensate for an incorrect 
reference cosmology i.e. the scaling relations given in Eq. (HJi, 
are sufficiently accurate for our simulated surveys. We calculate 
the distances assuming three reference cosmologies that differ 
only in terms of w,ef and compute the corresponding power spec- 
tra. 

As an example, the oscillatory parts of the dark matter power 
spectra at z = 3 in real space are shown in Fig.[TOl We observe 
that the scale of the oscillations is compressed and stretched for 
Wref = -1.2 and Wref = -0.8, respectively. 

We scale the one-dimensional power spectra according to the 
scaling relation, i.e. P(k) — > /tr*^^ P{k/Aiso), where the factor A^^^ 
is due to the scaled volume element. The fractional difference 
of these scaled power spectra with respect to that of the correct 
cosmology, i.e. w = -1, is shown in Fig. [TT] The data points 
for w - -1.2 (squares) and w - -0.8 (plus signs) are shifted 
by -1-0.25% and -0.40%, respectively, to center them on the zero 
line. These small shifts have a similar origin as in the light-cone 
versus snapshot comparison; we calculate the scaling factor at 



Constraining w with the (an)isotropic power spectrum 

The two-dimensional power spectrum provides the possibility to 
measure /ty and Aj_, instead of the combination Ai^^ only, as in the 
case of the one-dimensional power spectrum. Since both scaling 
factors depend on w in different ways, we expect that w is more 
accurately constrained by the two-dimensional power spectrum. 
In the case of a simple constant w, this is not the case: both the 
mean values and the errors are similar to those derived from the 
one-dimensional power spectra. For more complex models of w, 
the measurement of /l|| and A± is very helpful as we see from 
Fig. [12] where we applied the model w - wq + (I - a)wa- For 
display purposes, we performed a coordinate transformation of 
the parameters to the "pivot" system wo + (1 - a)wa - + 
(dp - a)wa where we chose Zp = 0.3 as discussed below. The 
left (right) panel indicates the constraints obtained from the dark 
matter light-cone power spectrum around z = 3 (z = 1 ) in real 
space. One sees that forz = 3 the contour lines are open towards 
negative Wa even for the two-dimensional case. Nevertheless, in 
both redshift cases the use of the anisotropic power spectrum 
tightens the constraints substantially. 
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Fig. 11. The fractional difference of the scaled real-space power spectra 
at redshift z = 3 is shown for the three different reference cosmologies. 
Note that the data points for w = -1.2 and w = -0.8 are shifted verti- 
cally by -^0.0025 and -0.0040, respectively, to enhance the visibility of 
the scatter. The lines show the intrinsic error in the power spectrum due 
to cosmic variance (dashed line) plus shot noise (solid line). 
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Fig. 12. This plot shows the 68% and 95% confidence levels for the 
redshift-dependent w model obtained from the one-dimensional (dashed 
lines) and two-dimensional (solid lines) dark matter power spectrum at 
redshift z = 3 (left) and z = 1 (right) in real space. 



Estimates for future surveys 

For a prediction of how well upcoming observations will en- 
able w to be measured, we analyze the mock galaxy catalogs 
introduced above, namely the weakly and strongly biased one 
million galaxy catalogs at redshift z - 3 and the one and two 
million galaxy catalogs at z = 1. In Fig. [13] we present typical 
fitting results of the scaling factors for each type of catalogs. For 
the solid lines, we used the aforementioned priors on the cosmo- 
logical parameters, whereas for the dashed lines we fixed these 
parameters. 

Since the error in the scaling factors is dominated by the er- 
rors in the power spectrum, improving the accuracy of the cos- 
mological parameters to higher than that of the assumed priors 
does not significantly reduce the error in the scaling factors. 
However, it tightens the constraints on w (see Fig. [14] for con- 
stant w): the uncertainty in Hq, in particular, degrades the mea- 
surement of w. 

For the strongly biased sample at z = 3 of one million 
galaxies in a volume of lOGpc^, the mear|3 uncertainties in the 
Hubble parameter {H(z) « A^) and angular diameter distance 
(Da(z) cc 1/^^) are 1.8% and 1.2% (68% c.l.), respectively. This 
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Fig. 13. Typical contour plots of the joint PDF of and Aj^ obtained 
from the mock catalogs at redshift z = 3 (left: weak (top) and strong 
(bottom) bias) and z = 1 (right: 1 million (top) and 2 million (bottom) 
galaxies) are shown. The contour lines correspond to 68% and 95% con- 
fidence level. The dashed line shows the results when the cosmological 
parameters are fixed at the values of the simulation. 



corresponds to an error of ~ 12% for a constant w. By keeping 
the cosmological parameters fixed, the uncertainty in w is low- 
ered to ~ 5%. 

The coiTesponding numbers for the two million galaxy mock 
catalog at z = 1 are 2.8% and 1 .6% for the Hubble parameter and 
angular diameter distance, respectively. In this case, we derive an 
accuracy of ~ 1 1 % and ~ 4% for a constant w with Planck priors 
and fixed cosmology, respectively. 




We used 10 samples to derive the mean values. 



Fig. 14. PDFs for constant w derived from the one million strongly 
biased galaxy sample at z = 3 (left) and from the two million galaxy at 
z = 1 (right). 

To translate the results for the scaling factors into constraints 
of the redshift-dependent w-model w = wo + (1 - a)wa, we com- 
bine the data of the strongly biased one million galaxy catalog 
at z = 3 and two million galaxy ca t alog at z - 1 with 192 
SN la observations (Davis et al. '2007^ IWood-Vasey et alJl2007l: 
Riess et al. 2007; Astier et al. 2006). In addition we include con- 
straints from CMB measurements by holding the ratio of the dis- 
tance to the last scattering surface and the sound horizon fixed: 
Da{z = zi,ss)/i = const. 

The confidence contours for wq and Wa are shown in the up- 
per panels of Fig. [15] In the lower panels, the confidence con- 
tours are presented in the pivot variable Wp = wq + (1 - ap)Wa, 
where we choose Op such that the error "ellipse" from the com- 
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bined datasets (solid line) is least tilted; we find Op — 1/(1 +0.3). 
One sees that the intersection angle of the BAO contours with 
respect to the error ellipse derived from SN data varies with red- 
shift (see also Fig.[T6]l. For z = 3, the orientation is very similar 
to the CMB contours, whereas for z = 1, it is more aligned with 
the SN ellipse. Measurements at both redshifts would therefore 
constrain the parameters of this model significantly more than 
two (or twice a good) measurements at a single redshift. 
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Fig. 15. The 68% and 95% confidence level for the parameters of the 
model w = wq + (1 - a)w„ obtained from SN (dotted), SN+CMB 
(dashed), and SN+CMB +BAO (solid). Results for the strongly biased 
one million galaxy sample are shown at z = 3 (left) and for the weakly 
biased two million galaxy sample at z = 1 (right). At the bottom, we 
show the same results but transformed to the pivot system Zp = 0.3. 

In Fig. [16] we show the constraints on w derived alone from 
BAO measurements at two different redshifts, namely z = 1 
(dashed lines) and z = 3 (dotted lines), and the combination of 
both measurements (solid lines). The redshift-dependent shape 
and orientation of the contour lines is clearly evident in the pivot 
system (right panel). 

For all confidence contours in Fig. [15] and Fig. [16] the above 
stated priors for the cosmological parameters Wm, cob, n^, and i/o 
were used. 

6. Discussion and Conclusions 

We have simulated large redshift surveys by "observing" mock 
galaxies on a light cone obtained from an N-body simulation 
with 512^' dark matter particles in a (1.5 Gpc/Zi)-' box. By fitting 
the apparent scale of the BAO in the power spectrum, we have 
derived constraints on the EOS parameter w. 

Our fitting method uses only the oscillatory part of the power 
spectrum and is therefore insensitive to changes in the power 
spectrum due to nonlinear evolution, redshift distortions, and 
scale-dependent galaxy bias. Fitting methods that use the overall 
shape of the power spectrum suffer in general from these compli- 
cated physical effects (unless they are accurately modeled) and 
tend to provide biased results (Smith et al. 2008). 

The drawback of fitting only the oscillatory part is that 
it is insensitive to information provided by the overall shape 
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Fig. 16. The 68% and 95% confidence level for the parameters of the 
model w = wo + (1 - a)w„ obtained from the BAO measurements in 
the strongly biased one million galaxy sample at z = 3 (dotted) and in 
the weakly biased two million galaxy sample at z = 1 (dashed). The 
solid lines show the combination of both BAO measurements. On the 
right, the contours are transformed to the pivot system of the combined 
dataset ap = 0.7. 



of the power spectrum. Determining the growth function 
would be interesting in particular for dark energy constraints 
(Anjendolaet al. 2005; Sapone & Amendola 2007). 

It is unclear how to extract the BAO from the power spectrum 
in the most appropriate and accurate way. Our simple smooth- 
ing method works fairly well. It is essentially parameter-free and 
produces unbiased and robust results. 

We analyzed the degeneracies in the EOS w with other cos- 
mological parameters. Given the expected accuracies in the cos- 
mological parameters involved, we found that more accurate 
measurements of //q would have the largest effect on lowering 
the uncertainties in w, especially for observations at lower red- 
shifts. 

Comparing the results of the light-cone power spectrum with 
those of the power spectrum of the snapshot at the correspond- 
ing redshift, we did not find evidence for substantial differences. 
For the surveys under consideration, light cone effects therefore 
do not play a role in determining w with our fitting method. This 
is in addition true, when we include redshift-dependent galaxy 
bias. Fitting methods, sensitive to the overall shape of the power 
spectrum, in particular to the growth function, must include ex- 
plicitly the light-cone effects to produce unbiased results. 

The PIZA technique ( Croft & Gaztanagal \l99% , which is 
used to reconstruct the linear density field and thereby the ampli- 
tude of the BAO, is effective if the density field is known to suffi- 
cient accuracy. In the case of substantial shot noise and unknown 
galaxy bias, more sophisticated reconstruction techniques are re- 
quired. 

To investigate the cosmological distortion due to an incorrect 
reference cosmology, we adopted three different reference cos- 
mologies with different values for Wi-ef = -0.8, -1.0, and -1.2. 
Employing the usual scaling relations (see Eq. (O), we found in 
all three cases similar fitting results within the margins of error. 
We conclude that even for this large survey volume the simple 
scaling of the power spectrum is fairly accurate. Nevertheless, 
for real data we propose to use an iterative scheme; after mea- 
suring w for an assumed reference cosmology, one would then 
repeat the analysis for the updated reference cosmology, using 
the measured value of w for Wref . 

By fitting the two-dimensional power spectrum, we can de- 
termine independently both scaling factors. For a constant w, the 
two-dimensional fitting does not improve the constraints on w 
significantly. For models of w that are not so tightly constrained, 
independent measurements of the Hubble parameter (//(z) cc A\\) 
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and the angular diameter distance (Da « 1 //t^) will however be 
very helpful. For the model w = wq + (1 - a)wa, this is especially 
true at lower redshifts. 

Our estimates for future surveys show that BAO measure- 
ments around redshifts z = 3 and z = 1 in combination with 
SN and CMB data tighten the constraints on dynamical dark en- 
ergy substantially and that these redshift surveys deliver comple- 
mentary results. For a constant w-model the BAO measurements 
from our two mock catalogs at z = 3 and z = 1 alone will con- 
strain Wo to an accuracy of ~ 12% and ~ 11%, respectively. 

To achieve more realistic predictions, the survey geometry, 
i.e. the window function, has to be taken into account. A more 
realistic galaxy bias scheme, tailored in particular for the tar- 
get galaxies, should also be developed. For such a study, large- 
volume simulations of far higher mass resolution are needed. 
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